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SUPERSONIC FLUTTER OF FLAT RECTANGULAR ORTHOTROPIC PANELS 
ELASTICALLY RESTRAINED AGAINST EDGE ROTATION 

By Larry L. Erickson 
Langley Research Center 

SUMMARY 

A theoretical solution is presented for the supersonic flutter characteristics of 
flat orthotropic panels subject to inplane loads and various boundary conditions. The 
solution is valid for panels having nondeflecting edges which are simply supported, 
clamped, or elastically restrained against rotation. 

Numerical flutter results of the analysis are tabulated in terms of general param- 
eters for simply supported and clamped panels and for three intermediate conditions of 
rotational edge restraint. From these parameters the dynamic pressure, frequency, and 
mode shape at flutter can be determined for large ranges of length-width ratio, stiffness 
ratio, and inplane stress. 


INTRODUCTION 

The flat rectangular panel is a basic structural element used for exterior skin sur- 
faces of vehicles which operate within the atmosphere. The dynamic instability, or 
flutter, of such panels when exposed to a supersonic flow has been the subject of numer- 
ous theoretical investigations. (See ref. 1 for a summary of panel flutter research.) It 
has been noted in previous studies that large differences exist between the flutter bound- 
aries of simply supported (ref. 2) and clamped (ref. 3) panels which are subjected to com- 
pressive stresses. These differences are especially prominent in the vicinity of the 
stress that produces buckling (ref. 4). Hence, there is a need to obtain flutter bound- 
aries for panels whose edges are supported in a manner intermediate to the simply sup- 
ported and clamped cases. To obtain these flutter boundaries, it is necessary to con- 
sider panel supports which offer resistance to edge rotation. 

The flutter solutions for flat rectangular panels with various degrees of rotational 
edge restraint are obtained in this report. The panel edges are assumed to be restrained 
from angular rotations by uniform restoring moments which are of equal strength on 
opposite edges. There is assumed to be no displacement of the panel edges in a direc- 
tion perpendicular to the plane of the panel. (Such displacements, if allowed to occur, 
can significantly affect the flutter boundaries as shown in ref. 5.) By accounting for the 



effect of rotational edge restraint, the results presented herein should be useful in 
obtaining more accurate flutter boundaries of compressively stressed panels, provided 
that the rotational restraint of the supporting structure can be estimated. 

SYMBOLS 

A parameter defined by equation (6) 

Aj constants appearing in equation (11); j = 1, 2, 3, 4 

a panel length in x- direction 

B parameter defined by equation (6) 

Bj,B 2 ,Bg,B 4 ,Bg,Bg coefficients defined by equation (A3) 
b panel width in y- direction 

Cq,Cj,C 2 coefficients defined by equation (9) 

D x ,Dy panel bending stiffness in x-direction and y- direction, respectively 
Dxy panel twisting stiffness 

Di,Di 2 ,D 2 panel stiffness coefficients defined by equation (8) 

F expression for the determinant appearing in equation (17) 

j,k integers 

N b 2 

k x nondimensional stress coefficient, — £ — 

7T 2 Dl 

Nyb 2 

ky nondimensional stress coefficient, — 

l aerodynamic pressure load per unit area given by equation (2) 

M Mach number 

mj roots of equation (12); j = 1, 2, 3, 4 


2 


N-v 


inplane loading in x- direction, positive in compression 


Ny 

q 

q x 

q y 

Sl,S 2 ,S3 


inplane loading in y- direction, positive in compression 

1 o 

dynamic pressure of airflow, -p a V 

£i 

rotational restraint coefficient on leading- and trailing- edge boundaries, 
a vx 


Dl 


be. 


rotational restraint coefficient on side- edge boundaries, — 

D 2 

coefficients defined by equation (16) 


time 


velocity of airflow 


w 


lateral deflection of panel 


X,Y functions of x only and y only, respectively, used in approximate solution 

to equation (1) 


x,y,z 


rectangular Cartesian coordinates (see fig. 1) 


a,6,e components of the roots mj (see eq. (13)) 

0 = /m 2 - 1 


e = ie 


spring constant of rotational springs supporting panel at x = 0 and x = a 


Mx> My 

Pa 


spring constant of rotational springs supporting panel at y = 0 and y = b 

dynamic pressure parameter, - 

Poisson's ratio in x- and y- directions, respectively 

mass density of airflow 
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Pm 


mass density of panel per unit area 


o> 


Subscripts: 

cr denotes critical or flutter value 

T denotes transition value (where frequency becomes zero and panel is on 

verge of buckling) 


panel frequency 
reference frequency, 


tt 4 Dj 

b 4 p m 


ANALYSIS 

The configuration analyzed consists of a flat rectangular orthotropic panel which 
is of length a and width b as shown in figure 1. The panel is subjected to uniform 
inplane force intensities N x and Ny which are considered positive in compression. 
The inplane shear intensity is taken to be zero. The panel is supported in such a man- 
ner that there is no lateral deflection along the edges. In addition, the edges are elasti- 
cally restrained against rotation by a uniform restoring moment which is proportional to 
the slope of the panel at the boundaries. The proportionality between restoring moment 
and slope is assumed to be equal on opposite edges. The supersonic flow at Mach num- 
ber M is over the top surface of the panel and is parallel to the X-axis. 


Differential Equation and Boundary Conditions 

The differential equation governing the lateral vibrations of a flat orthotropic panel 
is obtained from reference 6 by neglecting terms involving shear deformations and by 
adding a transverse inertia term 


D, 


1 — Mx My 9x^ 


^3 + 2 (° x y 


Myl^x 


8 4 w 


Dy 


9 4 W 


1 " MxM yl 9x 2 8y 2 1 - MxMy 9y4 


_ _ 9 2 \v d^W 9 2 W 

+ N x p- + Ny — =- + P m — p- - 

ax 2 ay 2 at 2 


l 


( 1 ) 


In this equation, D x and D y are the panel bending stiffnesses in the x- and 
y-directions, respectively; Djjy is the panel twisting stiffness; p x and p y are 
Poisson's ratios in the x- and y-directions, respectively; p m is the mass density per 


4 


unit area of the panel and l is the lateral pressure due to the airflow. The aero- 
dynamic loading is assumed to be given by two-dimensional static aerodynamics so that 


, _ -2q 9w 
0 ax 


( 2 ) 


where q = ^p a V 2 is the dynamic pressure and j3 = /m 2 - 1. Reference 7 has shown 
that for Mach numbers greater than about 1.5, the use of this simple expression to repre- 
sent the aerodynamic pressure load leads, in most cases, to the same flutter results as 
the use of more exact aerodynamics. 


The boundary conditions which the solution to equation (1) must satisfy are 
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9w 

9y 


Ay^=0 


and w = 0 


and w = 0 


and w = 0 


and w = 0 


at x = 0 I 

at x = a 
J 

"\ 

at y = 0 

> 

at y = b 


(3a) 


(3b) 


where 0 X and Ay are the spring constants (per unit length) of the rotational restraints 
acting at the boundaries. 


Solution of Differential Equation 
In general, a product solution of the form 

w(x,y,t) = x(|)Y^)e iwt (4) 

where go is the panel frequency, will not satisfy equation (1) since the term 

prevents the functions X and Y from separating. However, an approximate solution 
in the form of equation (4) can be obtained by use of a method given by Kantorovich 
(ref. 8). If Y(y/b) is assumed to be some function (as yet unspecified) which satisfies 
the boundary conditions given by equations (3b), then equation (1) can be reduced to an 
ordinary differential equation by the following procedure: 


A 
9 w 

9x3 9y 2 
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ii 


(1) In equation (1), w(x,y,t) and l are replaced with equations (4) and (2), 
respectively. 

(2) The resulting equation is multiplied by Y(y/b) and integrated with respect 
to y/b. 

This procedure yields 

X ""(i) + * 2Sx "(f) + *X’(|) - ir 4 Bx(|) = 0 (5) 

where a prime indicates differentiation with respect to x/ a and 



The coefficients Cq, Ci, and C 2 are due to the integrals involving Y(y/b) and are 
given by 
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where a prime indicates differentiation with respect to y/b. 


( 9 ) 


Thus, by selection of an appropriate function for Y(y/b), the problem is reduced 
to finding the exact solution of equation (5) which satisfies the following boundary 
conditions: 


where 


X(0) = X(l) = 0 
X"(0) - q x X'(0) = OS 
X"(l) + q x X’(l) = Oj 



is defined as the rotational restraint coefficient. 
The general solution to equation (5) is 


( 10 ) 


x(|) = Al e m la 


™ X 

+ A2e + A 


w. x 

3 e 3a + A^ 


m X 

, m 4a 


(ID 


where mj (j = 1, 2, 3, 4) satisfies the auxiliary equation 

m 4 + w 2 Am^ + Xm - ir 4 B = 0 (12) 

Equation (12) is of the same form as that solved by Hedgepeth (ref. 2) who assumed that 
the roots are 


mj = a + iS' 
m 2 - a- i6 
mg = -a + e 
m 4 - -ct- 


(13) 
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From the relations between the sum of the roots (taken one, two, three, and four at a 
time) and the coefficients of a polynomial equation, the following expressions are 
obtained: 


4a 

X_ 

4a 


Equation (14c) can be written as 



(14a) 



(14b) 

/ « \2 

wP 

M It}< 
+ 

1 

(14c) 

¥ S 2 a 2 - S 3 = 0 

(15) 


where 


q - J A 

S !-' T 

/ 2-\ 2 

4 b .IrA) 

7rB+ \^2~) 


S 2 =T 




(16) 


Thus, a can be determined from equation (15) for given values of the coefficients A, 

X, and B. The solution to equation (12) is then readily determined since 6 and e can 
now be calculated and hence the roots mj are known. 

Satisfaction of Boundary Conditions at the Leading and Trailing Edges 
Application of the boundary conditions, given by equations (10), to equation (11) 


gives 


m l( m i-q x ) m 2 (m 2 -q x ) 


, m l 


_ m 2 


m 3 ( m 3 " <hc) 
“3 


|mi (mj + q x )e mi m 2 (m 2 + q x )e‘" z m 3 (m 3 + m 4 (m 4 + q x )e“‘*| 


“2 


m 3 


m 4 (m 4 - q x ) 

e m 4 

, m 4| 


ft? 


■ 


a 4 


> (17) 


0 


y 


The condition for a nontrivial solution is obtained by equating the determinant of the 
square matrix to zero. This yields 
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/ O o\/ 9 o\f m 3 +m 4 m l +m 2 > \ 

F( m i, m 2 , m 3 ,m 4 ) = - m 3 ^m 2 ^ - m^e + e ) 

+ (m 4 ^ - m^^mj 2 - m 3 ^)( e m ^ +m ^ + e m l +m 3^ 

+ (V-V)K 2 --2 2 )(e mi+m4 + e m2+m3 ^ 

+ q x (m 2 - m 1 )(m 4 - m 3 ) [(mj + m 2 ) - (m 3 + m 4 ^e m3+m4 - e mi+m2 j 
+ q^n^ - m 3 )(m 4 - m^mj - m 2 ) + (m 3 - m 4 ^< 

+ q x (m 4 - m 1 )(m 3 - m 2 ) |m 2 - + (m 3 - m 4 )J(, 

“ t lx 2 ( m 4 " m 3)( m 2 " m 1 )(e m3+m4 + e 1 2 ) 

„ ( m 9 +m. m 1 +mo\ 

- q x 2 (m 4 - m 2 )(m 1 - m 3 ^e + e 3 ] 


e m 2 +m 4 _ e m 1 +m 3 j 

mj+m^ m 2+1113^ 

© — e I 


2 ( m l 

- <*x ( m 4 - m l)( m 3 - m 2^ e 


+m 4 m 2 +m 2 


+ e 


') =0 


(18) 


Replacement of the roots mj with the expressions given by equations (13) yields 

F(a,5,e) = [(e 2 + 6 2 ) 2 - 4a 2 (e 2 - 6 2 )]sinh e sin 6 - So^Sfcosh e cos 6 - cosh 2 a) 

+ 2q x [4a6e sinh 2a + e(e 2 + 6 2 - 4a 2 )cosh e sin 5 - s(e 2 + 6 2 + 4o 2 )sinh e cos e] 
+ q x 2 |2e6(cosh 2a - cosh e cos 6) + (e 2 - 5 2 - 4o?)sinh e sin = 0 (19) 


For the case where the panel is simply supported at x = (0,a), q x equals zero 
and equation (19) reduces to the solution obtained by Hedgepeth (ref. 2). For the case 
where the panel is clamped at x = (0,a), q x is infinite and equation (19) reduces to the 
solution obtained by Houbolt (ref. 3). 

Numerical calculations show that for values of A greater than about 4, equa- 
tion (14b) yields e 2 < 0. When this situation occurs it is convenient to define 

= ie so that equation (19) becomes 
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iia ai ■ ■ i 


FM,€) = [(s 2 - F 2 ) + 4a 2 (s 2 + e 2 )]sin F sin 6 - So^FSfcos F cos 6 - cosh 2a) 

+ 2q x [4oie 6 sinh 2a + e (fi 2 - F 2 - 4o 2 )cos F sin 6 - s(6 2 - F 2 + 4a 2 ) sin F cos sj 

+ q x 2 j2e 6(cosh 2a - cos F cos 6) - (fi 2 + F 2 + 4o 2 )sin F sin sj = 0 (20) 

The derivation of the panel mode shape associated with equation (19) is presented 
in the appendix along with a discussion of the modal characteristics. 


Flutter Condition and Procedure for Obtaining Numerical Results 

Equation (19) or (20) implicitly gives the panel frequency parameter B in terms of 
the dynamic pressure parameter x. When X reaches a certain critical value, desig- 
nated x cr , two of the panel frequencies coalesce. An increase in x above X cr causes 
B, and hence u> 2 , to become complex. Thus, one of the two square roots of o> 2 must 
possess a negative imaginary part which by equation (4) produces a divergent motion, 
termed flutter. The point where two frequencies coalesce is characterized by the con- 
dition d\/ 8B = 0 (see ref. 2) and corresponds to the values of X cr and B cr pre- 
sented in this report. 

For given values of the rotational restraint coefficient q x and the parameter A, 

.the following procedure was used to calculate the numerical flutter solutions presented 
herein: 

(1) A reasonable value was chosen for the parameter B and a trial value of X 
was selected. 

(2) The quantity a was calculated from equation (15). 

(3) The quantities 6 and e were then determined from equations (14a) and (14b). 

(4) Equation (19) or (20) was used to calculate F(o,6,e); a nonzero value indicated 
an incorrect choice of X. (From a programing aspect, it is more convenient to cal- 
culate F(a,5,e) from eq. (18), provided that a digital computer capable of working in 
complex arithmetic is available.) 

(5) The procedure was repeated with different values of X until the transcendental » 

equation F(o,6,e) = 0 was satisfied within a 0.01-percent error in X, thus giving a 

point on a X-B plot. 

By repeating the process for additional choices of B, the frequency loop (variation in X 
with b) was obtained for the given values of q x and A. The quantities X cr and 


i 
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B cr were then determined from the point on the frequency loop where 9x/ 9B = 0. A 
typical plot showing the first two frequency loops for the simply supported and clamped 
cases is shown in figure 2. 


Preflutter Solution 

When large negative values of the parameter A are considered, the exact solu- 
tion for X cr becomes difficult to determine numerically. However, as discussed in 
the next section, it is in this range of A that a simple algebraic solution can be used 
to obtain good approximations to both X cr and B cr . The algebraic solution is 
obtained by noting that for either simply supported (q x = 0) or clamped (q x = «°) panels, 
equation (19) is identically satisfied if 6 = 2kn and e = 2a, where k is an integer. 
This solution was first noticed by Movchan (ref. 9). With the use of these expressions 
for 6 and e , equations (14) yield the following simple algebraic solution for X and 
B which has been denoted as the preflutter solution 



For a given value of A, these two expressions determine the points where the fre- 
quency loops of the simply supported panel intersect those of the clamped panel. As 
shown in figure 2 the integer k indicates the loops which intersect. The lowest value 
of X is obtained by setting k = 1 ; thus, the preflutter equation for X with k = 1 
always yields X ^ x cr . Note that with k = 1, the preflutter solution is valid only for 
A 14. 


RESULTS AND DISCUSSION 

The primary results of the analysis are presented in table I for q x = 0, 2, 10, 40, 
and <». For each value of q x the values of X cr and B cr are tabulated for a wide 
range of the parameter A. The corresponding values of a cr are also given in table I 
so that mode shapes can be calculated from equation (A2). The variation of X cr and 
B cr with A is presented in figures 3 and 4, respectively, for A ^ -1. For Ii-1, 
the variation of X cr with A is presented in figures 5, 6, 7, 8, and 9 for qx = 0, 2, 10, 
40, and <*>, respectively. 
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In figure 3 the parameter A 


cr 1/S /fS, 


which does not contain the panel length a, 


is plotted against negative values of A for various values of q x . As -A becomes large 


(implying large length- width ratios if the inplane load N x is zero) the results become 


nearly constant. Because the parameter A, 




does not contain the panel length 


a, this behavior indicates that the flutter dynamic pressure has become independent of 
the panel length. At a value of A near - 100 the clamped (q x = °°) and the simply sup- 
ported (q x = 0) flutter boundaries have merged. Consequently, for A ^ -100, the 
amount of rotational restraint acting at the leading and trailing edges no longer has an 
effect on A cr . At a value of A near -300 the exact flutter boundary (for 0 ^ q x ^ °°) 
has merged with the preflutter solution (k = 1). Thus, regardless of the amount of rota- 
tional edge restraint acting at the leading and trailing edges, A cr can be obtained from 
the preflutter solution (eq. (21)) for any panel configuration having an A ? -300. Note 
that, for q x = 0, the preflutter solution is quite accurate for any negative A. 


In figure 4 the parameter B cr /(-A ) 3 (which also does not contain the panel 
length a) is plotted against negative values of A. This parameter also becomes con- 
stant as -A becomes large. Hence, for stress-free panels with large length-width 
ratios, the flutter frequency is independent of the panel length. 


Compressively Stressed Panels 

For panels which are subjected to compressive inplane force intensities in the flow 
direction, A can take on positive values. The flutter boundaries in this case are shown 
in figures 5 to 9 as plots of A cr 1//3 against A for values of A from -1 to 21. (The 
cube root of the flutter parameter is introduced merely for convenience in plotting.) Two 
boundaries are shown in each of these figures. The solid boundary results from the 
coalescence of the lowest two natural frequencies. The dashed boundary is obtained 
from the coalescence of the next two higher frequencies. For any value of A, the flutter 
value of A is determined by the curve that gives the smaller value of the ordinate. The 
preflutter solution (with k = 1) is included in figure 5 for comparison with the exact flut- 
ter boundary of a simply supported panel (q x = 0). The results show that the preflutter 
solution gives a reasonably good approximation to A cr for a panel with simply sup- 
ported leading and trailing edges whenever A § 3. 

The flutter boundaries in figures 5 to 9 indicate that the degree of rotational edge 
restraint can greatly affect the flutter behavior of compressively stressed panels. This 
influence of boundary conditions can be more clearly seen in figure 10 where A^ 1 / 3 
is plotted against the midplane stress ratio k X j/k x <p for a particular panel configura- 
tion and for several values of the rotational restraint coefficient q x . The panel is 
isotropic (Dj = D 12 = D 2 ) with a/b = 3, N y/ /N x = 1, and 8 x /8y = 1. The boundaries 
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shown were constructed from the numerical values of A cr and B cr presented in 
table I. The quantity is the transition value of k x ; that is, that value of k x 

which produces buckling (co = 0). Such a plot is often used in correlating experimental 
results when the ratio of the applied stress to the buckling stress can be estimated even 
though the magnitude of the stress is unknown. Experiments indicate that the transition 
point ^k x = k x x) * s the lowest point on the flutter boundary since a further increase in 
stress (panel buckled) usually has a stabilizing influence (ref. 4, fig. 14). For k x = 0, 
the edge support conditions affect the panel behavior in an expected way; that is, an 
increase in rotational restraint q x is accompanied by an increase in A cr . However, 
as k x /k X T increases, this behavior is completely reversed. For instance, when fixed 
values of k x ^/k x 'p are greater than about 0.4 the effect of stiffening the panel supports 
is to lower the dynamic pressure required to induce flutter. Further increases in 
k x /k X T result in a wide divergence of the boundaries. It should be noted however that 
the abscissa of figure 10 represents the inplane load as a percentage of the buckling load. 
This buckling load k x -p varies with qx as shown. Hence, for a given value of 
k x /k x T ^ 0) each panel configuration shown (q x = 0, 10, carries a different inplane 
load. 

For a given panel configuration, the variation of flutter frequency with stress is 
given in figure 11 where B cr is plotted against positive A for simply supported 
panels (fig. 11(a)) and clamped panels (fig. 11(b)). The solid lines correspond to the 
flutter frequencies associated with the solid flutter boundaries shown in figures 5 and 9. 
Superimposed on figures 11(a) and 11(b) are dashed lines which indicate panel natural 
frequencies (i.e., A = 0). For the range of A shown, the flutter frequencies lie 
between the two lowest natural frequencies. Note, however, that the two lowest frequen- 
cies do not always correspond to the first and second mode shapes. (See ref. 10 for a 
detailed discussion of this behavior.) Values of B cr for other values of q x can be 
obtained from table I; the curves shown in figures 11(a) and 11(b) are typical of the 
results for all values of q x . 

The theoretical results presented herein show that flutter is possible at zero 
dynamic pressure. It has been shown in reference 11 that, if damping is considered, 
flutter at zero dynamic pressure does not occur. However, the resulting theoretical 
flutter boundaries, which include the effect of damping, still do not agree well with 
experimental boundaries (ref. 7). The theoretical effect of damping then would be to 
smooth out, but not eliminate, the saw-tooth nature of the flutter boundaries in figures 5 
to 9. Thus, these flutter boundaries should not be considered reliable in the regions of 
positive A for which the theoretical values of A cr 


are near zero. 



Application of Results 


The results of the analysis show X cr to be a function only of A and the rota- 
tional restraint at the leading and trailing edges. The effect of the rotational restraint 
along the streamwise edges on the flutter solution is reflected in the coefficients C]/Cq 
and C 2 /C 0 which appear in A and B. Before the results tabulated in table I can be 
applied to a specific configuration, these coefficients must be determined. Values of 
these coefficients were obtained in reference 10 by the use of beam modes for the mode 
shape Y(y/b) appearing in equations (9). The beam modes considered satisfy the 
boundary conditions (3b). Thus, the resulting values of Cj^Cq and C 2 /C 0 are func- 
tions of the rotational edge restraint By along the boundaries parallel to the flow. The 
results for the first symmetrical mode are taken from reference 10 and presented in 
figures 12 and 13 where the variations in Cj^/Cq and C 2 /C 0 , respectively, with 

_ b0£ are s h own> Only the solution for the first mode is shown here since it usually 
' D 2 

yields the most critical flutter condition. (Ref. 10 presents solutions for the first four 
modes.) 

Once the ratio C^yCg is known, the parameter A can be determined from equa- 
tion (6a) for a given panel configuration and stress condition. Then, for a specified 
amount of rotational edge restraint along the boundaries perpendicular to the flow (given 
by qx), the flutter value of X and B can be determined from table I or from the fig- 
ures showing X cr and B cr as a function of A and q x . The flutter frequency can 
then be determined from equation (6c). It is important to make this calculation for the 
frequency because a nonpositive value of indicates that the panel is buckled, 

and the flutter boundaries presented herein are no longer valid. 

Throughout this section it has been implied that the spring constants 0 X and By- 
are known. These constants will, of course, depend on the configuration of the structure 
which supports the panel. The values can be obtained either from a stiffness analysis of 
the supporting structure or from experimental vibration frequencies of the supported 
panels. 

It should be remembered that the solutions presented herein are approximate 
unless the side edges (y = 0,b) are simply supported, in which case they are exact. How- 
ever, the vibration and buckling results of reference 10, which were obtained by using 
the same approximation as used herein, were in good agreement with exact solutions. 
Therefore, except in cases where consideration of damping and exact aerodynamics is 
important, the flutter results presented in this report are expected to be reasonably 
accurate. 
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CONCLUSIONS 


The theoretical solution is obtained for the flutter behavior of a flat orthotropic 
rectangular panel supported in such a manner that the rotations of the panel edges are 
elastically restrained. The solution is in terms of general parameters which account 
for dynamic pressure, Mach number, panel frequency, length-width ratio, inplane stress, 
panel stiffnesses, and edge support conditions. Numerical flutter values of these param- 
eters are tabulated for simply supported and clamped panels and for three intermediate 
conditions of rotational edge restraint. From the results of the analysis, the following 
conclusions can be drawn: 

1. The flutter boundary for compressively stressed panels is very sensitive to the 
panel support conditions. 

2. When the length-width ratio is large the flutter dynamic pressure and frequency 
for stress-free panels are insensitive to the support conditions at the leading and 
trailing edges and the panel length. 

3. The simple algebraic preflutter expression can be used to calculate the flutter 
boundary for a wide range of parameters with good accuracy. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 23, 1966. 
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APPENDIX 


MODAL BEHAVIOR 

The panel mode shape can be obtained from equation (11) once the coefficients Aj 
are known. These coefficients can be obtained, within an arbitrary constant C, from 
equation (17) 
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Substitution of these expressions for Aj into equation (11) and replacing the mj with 
the expressions given by equations (13) yields, after considerable manipulation, 



4iC 


x 

ot— 


(Bi + q x B 4 ) e &sin 6 S + ( b 2 + q x B 5 ) e 


\ 


- a 5 

a cosh e | 



+ 


( b 3 + q x B e) 



sinh e ^ 

d 


(A2) 
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where 


®1 

b 2 

b 3 

b 4 

b 5 

b 6 


2a!e (e^cos 6 - e -Q! cosh e) - (e^ + 6^)e -Q! sinh 

2o'ee cx sin 6 + 2a!6e -Q! sinh e 

(e^ + 5^)e a sin 6 + 2a5 (e^cos 6 - e -Q! cosh e) 

e(e Q! cos 6 - e -a! cosh e) - 2o!e -Q! sinh e 

ee^sin 6 - 6e -Q! sinh e 

2o!e Q! sin 6 - 6(e Q! cos 6 - e -Q! cosh e) 


,1 




(A3) 


Equation (A2) satisfies the first three boundary conditions (eqs. (10)) identically and 
satisfies the fourth only if the transcendental equation F(af,6,e) = 0 is satisfied. 

The theoretical mode shape at flutter is strongly influenced by compressive values 
of the loading N x . This effect is indicated by the following table which shows the num- 
ber of nodes and the region where the maximum amplitude occurs for simply supported 
and clamped panels. (Increasing A reflects increasing compressive values of N x .) 


A 

Maximum amplitude 
occurs toward — 

Number of nodes 


Simply supported 

Clamped 

Simply supported 

Clamped 

-10 

Trailing edge 

Trailing edge 

1 

1 

0 

Trailing edge 

Trailing edge 

1 

1 

3 

Trailing edge 

Trailing edge 

1 

1 

6 

Leading edge 

Trailing edge 

1 

1 

9 

Leading edge 

Trailing edge 

1 

2 

12 

Trailing edge 

Leading edge 

2 

1 


Correlation of the information presented in this table with figures 11(a) and 11(b) 
shows that the point of maximum amplitude moves from the rear of the panel to the 
front and back again as A passes the points where natural frequency lines cross. The 
flutter mode shapes at A = 12 are shown in figure 14 for three different values of q x . 
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TABLE I.- FLUTTER SOLUTIONS FOR FLAT RECTANGULAR PANELS WITH 


VARIOUS DEGREES OF ROTATIONAL EDGE RESTRAINT 


A 


j? 

ii 

o 



qx = 2 



ii 

O 

A 

\ C r 

B cr 


*-cr 

B cr 


X C r 

B cr 

-300 

92,230 

8,584 

21.96 

92,240 

8,594 

21.91 

92,300 

8,606 

-200 

51,390 

4,060 

17.89 

51,410 

i 4,063 

17.90 

51,470 

4,078 

-100 

19,410 

1,200 

12.68 

19,430 

1,202 

12.70 

19,500 

1,212 

-70 

11,980 

667.5 

10.64 

12,000 

670.0 

10.66 

12,080 

680 . ( 



7,720 

5,830 

4,119 

2,608 

1,330 

794.6 

697.1 

603.1 

512.6 

426.0 

343.3 

264.9 

190.9 

121.8 

57.98 

0 

51.27 

94.54 

127.5 

145.5 

137.3 
99.85 

51.01 

0 

49.67 

96.08 

137.7 

172.7 

199.0 

212.8 

208.2 
183.2 


395.0 

285.0 

190.5 

113.5 

53.50 

30.05 

25.75 
21.80 

18.00 
14.30 

10.75 
7.500 

4.375 

1.375 
- 1.412 
-4 

- 6.400 

- 8.625 

- 10.69 

- 12.75 

- 16.25 

- 22.50 

- 29.37 

-36 

- 42.40 

- 48.55 

- 54.50 

- 60.39 

- 66.25 

- 72.60 

- 81.00 

- 92.60 


9.082 

8.164 
7.177 
6.014 
4.582 
3.669 
3.467 
3.242 
3.004 
2.752 
2.480 

2.165 
1.810 
1.387 

.8262 

0 

.7211 

1.047 
1.171 
1.163 

1.048 
.8293 
.4818 

0 

.4682 

.7857 

.9802 

1.093 

1.140 

1.132 

1.081 

.9919 


7,747 

5,858 

4,149 

2,640 

1,364 

830.5 

733.3 

639.6 

549.5 

463.1 

380.8 

302.6 

228.9 

160.0 

96.21 

38.16 

13.48 
57.70 
92.94 

116.6 

123.8 

108.6 

74.48 
32.25 
11.61 
53.91 
92.63 

126.1 

152.6 

169.8 

175.2 

166.3 


397.5 

287.5 

193.5 

116.5 

56.50 

32.50 
28.25 
24.21 
20.40 
16.55 
13*00 

9.625 

6.375 

3.375 
.4500 

- 2.250 

- 4.800 

- 7.212 

- 9.500 

- 11.75 

- 14.62 

- 18.87 

- 25.00 

- 31.60 

- 38.12 

- 44.55 

- 50.86 

- 57.10 

- 63.49 

- 70.20 

- 77.75 

- 87.00 


9.099 

8.183 

7.179 

6.014 

4.574 

3.686 

3.482 

3.264 

3.027 

2.788 

2.520 

2.224 

1.895 

1.508 

1.049 

.4702 

.1659 

.6184 

.8490 

.9288 

.9031 

.7715 

.5506 

.2484 

.0891 

.3915 

.6185 

.7700 

.8594 

.8941 

.8803 

.8234 
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TABLE L- FLUTTER SOLUTIONS FOR FLAT RECTANGULAR PANELS WITH 
VARIOUS DEGREES OF ROTATIONAL EDGE RESTRAINT - Concluded 




o 

II 

£ 



q x =« 


A 

^cr 

®cr 

“ cr 

^cr 

®cr 

«cr 

-300 

92,400 

8,638 

21.86 

92,630 

8,684 

21.95 

-200 

51,590 

4,100 

17.89 

51,790 

4,142 

17.89 

-100 

19,620 

1,232 

12.69 

19,790 

1,257 

12.71 

-70 

12,200 

695.0 

10.69 

12,350 

717.5 

10.67 

-50 

7,946 

422.5 

9.082 

8,085 

437.5 

9.128 

-40 

6,058 

310.0 

8.190 

6,189 

323.7 

8.197 

-30 

4,348 

215.0 

7.169 

4,470 

226.0 

7.226 

-20 

2,839 

134.5 

6.054 

2,949 

144.0 

6.097 

-10 

1,560 

71.00 

4.679 

1,655 

78.50 

4.723 

-5 

1,022 

45.50 

3.828 

1,107 

51.50 

3.894 

-4 

923.4 

41.00 

3.633 

1,006 

46.50 

3.712 

-3 

828.2 

36.50 

3.435 

908.6 

41.69 

3.521 

-2 

736.5 

32.20 

3.227 

814.5 

37.20 

3.316 

-1 

648.3 

27.90 

3.013 

723.8 

32.49 

3.115 

0 

563.8 

23.75 

2.788 

636.6 

28.25 

2.891 

1 

483.0 

20.02 

2.536 

553.0 

24.00 

2.660 

2 

406.3 

16.25 

2.276 

473.3 

19.75 

2.420 

3 

333.9 

12.50 

2.002 

397.6 

15.75 

2.160 

4 

265.8 

9.000 

1.700 

326.1 

12.00 

1.877 

5 

202.4 

5.625 

1.374 

258.9 

8.250 

1.578 

6 

144.1 

2.250 

1.029 

196.5 

4.750 

1.253 

7 

91.17 

-.9000 

.6670 

138.9 

1.250 

.9172 

8 

44.08 

-4.062 

.3209 

86.74 

-2.217 

.5804 

9 

3.373 

-7.170 

.0236 

40.27 

-5.614 

! .2654 

10 

30.27 

-10.30 

.1995 

0 

-9 

0 

11 

56.07 

-13.60 

.3453 

33.48 

-12.50 

.2005 

12 

73.14 

-17.00 

.4203 

59.51 

-16.00 

.3348 

13 

80.73 

-21.00 

.4384 

77.40 

-19.85 

.4102 

14 

78.58 

-25.69 j 

.4082 

86.54 

-24.10 

.4347 

15 

67.60 

-31.20 

.3398 

86.74 

-29.00 

.4169 

16 

49.96 

-37.45 

.2447 

78.60 

-34.70 

.3651 

17 

28.40 

-44.30 

.1359 

63.69 

-41.20 

.2882 

18 

5.411 

-51.56 

.0253 

44.20 

-48.32 

.1954 

19 

17.04 

-59.16 

.0774 

22.34 

-55.97 

.0966 

20 

37.48 

-67.05 

.1651 

0 

-64 

0 

21 

54.80 

-75.25 

.2339 

21.35 

-72.36 

.0878 
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Figure 1.- Panel geometry, coordinate system, and flow direction. 



Figure 2.- First two frequency loops for clamped and simply supported panels. A = 0. 
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Flutter frequency 

Natural frequencies 

(b) Clamped leading and trailing edges. 
Figure 11.- Concluded. 
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